Biexciton stability in carbon nanotubes 
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We have applied the quantum Monte Carlo method and tight-binding modelling to calculate the 
binding energy of biexcitons in semiconductor carbon nanotubes for a wide range of diameters and 
chiralities. For typical nanotube diameters we find that biexciton binding energies are much larger 
than previously predicted from variational methods, which easily brings the biexciton binding energy 
above the room temperature threshold. 
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Size-dependent optical excitations in nano-structures 
are at the heart of fundamental studies as well as conceiv- 
able applications Carbon nanotubes (CNTs) make no 
exception, showing very sensitive electronic and optical 
properties to the atomic structure and spanning a wide 
range of wavelengths . The stability of the excitonic 
states, neutral or charged optically excited electron-hole 
complexes, is determined by their binding energy with re- 
spect to thermal fluctuations. In quasi- ID systems, the 
binding energy can be much larger than in systems of 
higher dimensionality. In inorganic semiconductor het- 
erostructures, for instance, the exciton binding energy 
Ef is substantially larger than [3, [3, [1] the binding en- 
ergy 4 Ry* in a strictly 2D system [y] , where Ry* is the 
Rydberg energy of the host material. Analogously, the 
biexciton binding energy E£ x of two electron-hole pairs, 
optically excited in a two-photonprocess, is not limited 
to its 2D value of 0.77 Ry* 

Since CNTs are quasi-lD systems, obtained by rolling 
up a graphene sheet llj , they are characterized by rather 
large binding energies 12l . H^ . analogously to conju- 



gated polymers [14|, [lfj. On the other hand, in CNTs 



one expects strongly diameter-dependent binding ener- 
gies. Indeed, in addition to the increase of the ratio 
/ Ry* with decreasing diameter, due to the transi- 
tion from a quasi-2D system to a quasi-lD system, also 
electron and hole effective masses, which determine Ry* , 
change with the CNT diameter. Furthermore, due to 
the involved energy scales, in CNTs not only excitons 
but also biexcitons might be stable against thermal fluc- 
tuations at room temperature; also, the energy separa- 
tion can be larger then the linewidth, and optical de- 
tection of biexcitons should be possible. Contrary to 
inorganic semiconductors and semiconductor nanostruc- 
tures, where biexcitons have received considerable inter- 
est, there is only very little work devoted to biexcitons in 
CNTs. Following the pioneering work of Ando [l6j], the 
exciton binding energy has been calculated for several 
nanotubes, with different diameter and chirality, both 
within ab initio approaches 
methods 0, S El, 0, 13 



biexciton binding energy, which is presently not acces- 
sible to the more accurate first principles methods, has 
been computed only via an approximate variational ap- 
proach [25j . 

In this paper we use the quantum Monte Carlo (QMC) 
method to calculate the exact (dimensionless) binding 
energy E£ x /Ry* of biexcitons confined to the surface 
of a cylinder of diameter D. We find that the biexci- 
ton is much more stable and its binding energy much 
larger than estimated from variational methods, partic- 
ularly for intermediate to large D/a* B . Assuming homo- 
geneous dielectric screening and tight-binding estimates 
for the CNT effective masses, we also estimate E^* for 
several families of CNTs. We find that for realistic val- 
ues of the dielectric constant, Ef x can be comparable or 
larger than fcsT room even for the larger CNT diameters. 

Let us consider electrons and holes confined to the sur- 
face of an infinitely long cylinder of diameter D, as indi- 
cated in Fig.[TJ It is convenient to use dimensionless exci- 
ton units, where H = 1, and in which masses are measured 
in units of the reduced electron-hole mass fi, distances in 
units of the effective Bohr radius a* B = (e/fj,)aB, and en- 
ergies in units of the effective Rydberg Ry* — e 2 /(2ea* B ). 
With the approximation of equal electron and hole ef- 
fective masses lj)|, 21, 27 1, the biexciton Hamiltonian in 
dimensionless units reads 



H = 



17, [1_8| and semiempirical 
24j . On the contrary, the 
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Here, r\,r 2 (r a ,ri,) are the positions of the two electrons 
(holes), and m is the distance between particles i and 
j. We choose electrons (holes) with opposite spin ori- 
entations and consider the optically active spin-singlet 
biexciton ground state. 

Variational QMC. — In the following we sketch our nu- 
merical approach. In addition to exact QMC calcula- 
tions, to be discussed below, we have performed vari- 
ational QMC (VQMC) calculations. We have exploited 
the (unnormalized) Hylleraas-Ore trial wavefunction [28| 
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FIG. 1: Biexciton complex on a cylindrical surface, x is the 
circumferential and z the longitudinal direction of the tube 
with diameter D = 2R. Particles 1,2 (electrons) and a,b 
(holes) are confined to the surface and form a spin singlet 
state. 



in a slightly modified version 



*T = e -^- + s "+s-+^) COsh[^( Sla - 8lb + S 2b - S 2a )\ . 

. (2) 

Here, the s^-'s are relative distances scaled by variational 
parameters. As we employ cylindrical coordinates, the 
expression for Sij reads 




D sin (xij /D) 



q 



(3) 



where x and z are oriented along the circumference 
and the symmetry axis of the cylinder, respectively (see 
Fig. [IJ. Variational parameters q and k allow for dif- 
ferent scaling for both directions [HI]. P is an addi- 
tional variational parameter determining the strength of 
the coupling between the two excitonic complexes of a 
biexciton, where (3 = 1 corresponds to two separate 
excitons, and for (3 — there is equal binding within 
all electron-hole pairs. The variational parameters q, k 
and (3 have to be determined such that the total energy 
E T = | If | *)/(*|*) becomes minimized. 

Quite generally, after separation of the biexciton 
center-of-mass motion, the calculation of Et involves 
six- fold integrals, which constitutes a formidable com- 
putational task. The calculation of Et is performed by 
the VQMC approach j^, whose main elements can be 
summarized as follows: since the trial wavefunction ([2|) 
of the optically active spin-singlet biexciton groundstate 
is always positive (thereby avoiding the fermionic sign 
problem), it can be represented by an ensemble of "walk- 
ers" , each one characterized by the particle positions r\ , 
r 2, fa, fb- Starting from a suitable initial configura- 
tion, one generates a Markov chain for the walkers where 
the probability for a specific configuration is given by 
^ T (ri,r 2 ,r a ,rb)- Upon sampling of the "local energy" 
El = H^t/^t one then obtains the energy Et asso- 
ciated to the trial wavefunction [29!]. Let us denote the 
ensemble of walkers with p{r\,r2,r a ,rb,t), where t is a 
fictitious time. The Fokker-Planck equation, which in 



our dimensionless units reads 
dp 1 



Of 



i=\.2,a,b 



V, (V, - Fi)p, 



(4) 



in time-discretized form defines a scheme to proceed from 
a configuration pit") to p(t + 5t) according to the drift and 
diffusion process given on the right-hand side [3(J. Thus, 
the VQMC simulation consists of the three main steps 
of (i) initialization of the ensemble of walkers, (ii) drift 
and diffusion of all particles in each walker according to 
Eq. 0, and (iii) sampling of the local energy El once 
the stationary distribution is reached. 

Guide Function QMC. — A slight variant of the VQMC 
approach allows for the exact solution of the Schrodinger 
equation. Let p = 'F'Ft denote a function composed of 
the exact wavefunction 'F and the guide function 'Fy. 
The Fokker-Planck-type equation 



dp 
dt 



V i (V i -F i )p-[EL-E]p (5) 



i=l,2,a,b 



in time-discretized form again defines a scheme that can 
be solved by means of Monte-Carlo sampling. It can 
be easily proven that under stationary conditions 
Eq. ((5l) reduces to the exact Schrodiger equation. Thus, 
once the invariant p is obtained the exact wavefunction 
is at hand. We represent p by an ensemble of walkers, 
and account on the right-hand side of Eq. ([5]) for the 
first term through drift and diffusion and for the second 
term through a Monte-Carlo branching with probability 
p = cxp[— (El — E)8t\. Depending on the value of p, the 
walker dies, survives, or gives birth to other walkers [29l |. 
In the simulation the energy E is chosen such that the 
total number of walkers remains approximately constant 
and a constant distribution p is reached. This invariant 
distribution and E then determine the biexciton wave- 
function and energy, respectively. Therefore, the main 
steps of this so-called guide function QMC approach are 
(i) initialization of the ensemble of walkers, (ii) drift and 
diffusion of all particles in each walker, (iii) branching of 
the walkers, and (iv) sampling of the wavefunction once 
the stationary distribution is reached. 

Technically, one needs to choose St sufficiently small 
to allow for the separate drift-diffusion and branching 
steps accounting for the two different terms on the right- 
hand side of Eq. ((5|) 3l| ; 'Ft has to be chosen such that 
the local energy El in Eq. ((5|) remains finite when two 
particles in a walker approach each other. While this 
is guaranteed for the true wavefunction, 'Ft is usually 
taken as a Jastro w-ty pe wavefunction with the correct 
cusp condition [29|, [34] . In practice we use ^t = with 
q = k = j, [3 = Q. We finally emphasize that both the 
VQMC and the exact QMC simulations can be applied in 
a straightforward manner to excitons, in which case the 
trial wavefunction is of the form 'F^ri, r a ) = exp[— s la ] 
and ^t = exp[— 2ri a ], respectively. 
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FIG. 2: Biexciton binding energies as a function of CNT diam- 
eter, calculated with exact guide function QMC (black dots) 
and VQMC method (grey triangles). Statistical error bars of 
both methods are shown. The solid, black line is the result of 
using our fitting functions in Eq. iff}. Inset: total energy of 
exciton and biexciton. Note that here E x = —E*. 



Results. — In the inset of Fig. [5] we show the dimension- 
less total energies E x and E xx with respect to the band 
gap (Eg — 0), calculated exactly by the guide function 
QMC method. As can be seen, the exciton energy E x 
shows the correct behaviour at infinite diameter, that is 
the 2D limit, where it converges to — 4i?y* 6], while it is 
strongly red-shifted as the diameter is decreased due to 
the larger binding of the electron-hole pair, showing the 
transition from a quasi-2D to a quasi-lD system. Analo- 
gously, the biexciton energy E xx red-shifts as a result of 
the increased interaction of the two electron-hole pairs. 

In order to investigate the stability of the biexciton 
complex with respect to the formation of separated ex- 
citons, we show in Fig. [2] the exact biexciton binding 
energy Ef x for arbitrary dimensionless diameter D/a* B , 
along with the VQMC results and the fitting function 
of Ref. l25|. The biexciton results to be stable (Ef x > 
0) at any diameter, and it shows the expected limit- 
ing behaviour at infinite diameter, that is the 2D limit 
0.77 Ry* The binding energy increases with decreas- 
ing diameter, showing again the transition from a quasi- 
2D to a quasi-lD system. As can be noted, both varia- 
tional results severely underestimate the binding energy 
in the whole range of diameters except for very small 
values. Moreover, they do not show the correct 2D 
limiting behavior. Such shortcoming of Hylleraas-Ore- 
type wavefunctions is in agreement with corresponding 
calculations for two-dimensional quantum wells 0, l33j ]. 
To give a rough estimate, for typical CNT diameters of 
D = 0.8-^1.2 nm and using the dielectric constant e = 3.5 
given in Ref. the excitonic units arc a* B = 2.5 -j- 5.5 
0.08 eV. This brings E xx into the 



nm and Ry* = 0.04 
0.06 -j- 0.12 eV range, which is 1.5 -j- 2.5 larger than the 
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FIG. 3: Exciton (a) and biexciton (b) binding energies as 
a function of CNT diameter and chirality. Inset: biexciton 
energies (—E xx ) in absolute units. The labels indicate the 
(2n+m) families, where (n,m) are the chiral indices of the 
NT. 



variational results reported in Fig. [2j 

In order to calculate exciton (x) and biexciton (xx) ab- 
solute energies explicitly, we assume that excitonic effects 
(x or xx binding) can be decoupled from band structure 
effects. We can therefore write the energies as 



EX = fx(D/a B )Ry*(D), 



(6) 



where \ G { x i xx } an d f x (x) is the exact dimension- 
less excitonic or biexcitonic energy shown in the inset 
of Fig.0 We provide below a fitting function for f x (x), 
which allows us to calculate absolute binding energies for 
an arbitrary diameter, 

f x (x) = (a x x^ + b x x~ 2 + c x x~ 3 ) exp(-d x x) + f x D , 

(7) 

where f x D is the correct 2D limit, and a x ,b Xl c x , d x are 
the fitting parameters summarized in table I. The quality 
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TABLE I: Fitting parameters in Eq. © for exciton (X) and 
biexciton (XX) energies in dimensionless exciton units. 





a x 


K 




d x 


fX 


X 


-2.62 


0.3024 


-0.01504 


2.795 


-4.00 


XX 


-5.08 


0.56 


-0.02832 


2.345 


-8.77 



of the fit is proven by calculating the biexciton binding 
energy with 2/ x — / xx , see solid black line in Fig. [21 

We calculate the Rydberg energy Ry* from the tight- 
binding model of Ref. I2R which provides explicit fitting 
functions [Eq. (2) of Ref. [2(| for the electron and hole 
effective masses of semiconducting CNTs of arbitrary 
chirality and diameter. For the dielectric screenig en- 
tering Ry* , we adopt the simplest screening model, in 
which all Coulomb interactions are reduced to an effec- 
tive static dielectric constant e, following current litera- 
ture [13L Il9l . |20 . 1 2 11 ] . This approximation, which proved 
to be successful in comparison with both experiments 34 1 
and ab-initio calculations [23TJ , is particularly suitable 
for embedding media with large dielectric constant, i.e. 
e > 3, where the dielectric response is dominated by that 
of the medium rather then by the CNT polarizability [13] • 

Figure [3] shows Kataura-like plots for both exciton (a) 
and biexciton (b) binding energies in absolute units, in 
the 0.5 -j- 1.5-nm diameter range, calculated from Eq. ([6]) 
and for different value of e. As expected, both x and xx 
energies decrease with increasing tube diameter, as fol- 
lows from the behaviour of f x (x). The chirality effects, 
entering through the effective masses, introduce a modu- 
lation in the binding energy dependence on the diameter, 
significantly spreading out the binding energies for the 
range of diameters considered. As shown in Fig. [3jb), 
however, the biexciton binding energy is predicted to be 
above the kBT rooal threshold (26 meV) even for the largest 
CNTs. 

In summary, we have performed QMC calculations of 
the singlet optically active biexciton binding energy for 
CNTs of arbitrary diameter and chirality. The biexciton 
has been found to be always stable at room temperature 
in typical dielectric environments. We have also devel- 
oped a scheme which allows us to calculate the exact 
exciton and biexciton binding energies through a sim- 
ple fitting function approach, which includes the strong 
correlation effects exactly and the band-structure effects 
within a tight-binding approach. 
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